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Abstract 



We present the analytical a^ correction to the Z° decay rate into hadrons. We 

calculate this correction up to (and including) terms of the order (w^/m^ ) 3 in the 

large top quark mass expansion. We rely on the technique of the large mass expansion 

,— i- J of individual Feynman diagrams and treat its application in detail. We convert the 

C^ ' obtained results of six flavour QCD to the results in the effective theory with five 

active flavours, checking the decoupling relation of the QCD coupling constant. We 

Qh! also derive the large charm quark mass expansion of the semihadronic r lepton decay 

O !* rate in the a^ approximation. 
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1 Introduction 

Precision measurements of the Z° decay rate into hadrons at LEP [|TJ provide precise means 
to extract the QCD coupling constant from experiment. This is a very clean process from 
a theoretical point of view since its calculation can be reduced to the calculation of the 
Z boson propagator within the standard model. The status of electroweak corrections to 
Z° decay can be found in ref . . Now the calculational techniques of Feynman diagrams 
have advanced so far that the calculation of the af order (=4 loop approximation of the Z 
boson propagator) is feasible. The af approximation to the Z° decay rate into hadrons is 
important for an accurate determination of the QCD coupling constant a s , or equivalently 
the fundamental scale of QCD, Aq C d- 

The hadronic Z° decay rate is a sum of vector and axial vector contributions of which the 
vector contribution is known to order a 3 from the calculation of <7t ((e + e~ — > 7 — > hadrons) 
@. This calculation was performed in the approximation of effective QCD with five mass- 
less quarks which involved the calculation of only massless diagrams. The correctness of 
this calculation is strongly supported by jp where the non-trivial connection between the 
result || and the a 3 approximation H to deep inelastic sum rules was established. 

The calculation of the axial vector part of the hadronic Z° decay rate is more involved 
than that of the vector part. This is because the heavy quark does not decouple in the 
axial vector part and one cannot avoid to calculate massive diagrams, even in the leading 
order of the large mass expansion. The axial vector part was calculated to order a 2 s in || 
and confirmed in |7[ where the operator product expansion technique was used to sum up 
the massive logarithmic terms. The Z° decay into 3 gluons in order a 3 has been calculated 
in ||. The af correction to the axial vector part of the hadronic Z° decay rate in the 
leading order of the large top mass expansion was presented in || [1(J. In this paper we 
elaborate the details of the calculation || and extend the large mass expansion of both the 
vector and axial vector parts to the order (m| / 'm^ op ) 3 . This calculation allows us also to 
check the decoupling mechanism at the next-next-to-leading order. 

Another prominent process (beside Z boson decay) for the extraction of a s from ex- 
periment is semihadronic r lepton decay. In the last section of our article we convert 
the obtained result for T(Z° — > hadrons) to the large charm quark mass expansion for 
r(r~ — > u T + hadrons) in the order a 3 ,. 

2 Preliminaries 

For the Z° decay rate into hadrons, the quantity to be determined is the squared matrix 
element summed over all final hadronic states. One can express this quantity as the 



imaginary part of a current correlator in the standard way 

Y < 0\J"\h >< h\J u \0 > = 2ImW\ 
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Here J M = ^-^ — Y%=\ ipij 1 * (dv ~ 9a1 5 ) & * s the neutral weak quark current coupled to 

the Z° boson in the Lagrangian of the Standard Model, where we use the notations as 
given in [JT1J g % v = t 3L (i) - 2% sin 2 W and g\ = t 3L (i). 
The hadronic Z° decay width is expressed as 
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with the indicated decomposition into vector and axial vector parts imposed by the struc- 
ture of the neutral current. We will calculate ImUi in the order g 2 a^. It is a calculation 
within perturbative QCD except for two weak current vertex insertions (i.e. the weak 
current is considered as an external current for QCD). 

Throughout this paper we use dimensional regularization |12| in D = 4 — 2e space-time 
dimensions and the standard modification of the minimal subtraction scheme |13| , the MS 
scheme [Q . For the treatment of the 75 matrix in dimensional regularization we use the 
technique described in [15] which is based on the original definition of 75 in [12]. We work 



in the approximation of 5 massless quark flavours and the top quark mass large compared 
to the Z° mass. We should stress that the top quark does not decouple |^6j from the axial 
vector part due to diagrams of the axial anomaly type. 

It is convenient to split the vector and axial vector contribution in non-singlet and 
singlet parts 
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The non-singlet parts come from Feynman diagrams where both weak current vertices are 
located in one fermion loop. The singlet contributions come from diagrams where each weak 
current vertex is located in a separate fermion loop. The massive non-singlet diagrams are 
presented in figure 1, the singlet diagrams are presented in figures 2,3. 

The o? s approximation for the vector part in effective QCD with 5 active massless quark 
flavours in the MS scheme was calculated in || (in the leading order of the large top quark 
mass expansion). This calculation used the fact that the top quark decouples for the vector 
part in the leading order of the large quark mass expansion. Therefore within effective 5 
flavour QCD this calculation involved only massless diagrams and the result reads 
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Here af\mz) is the coupling constant in 



effective QCD with 5 active flavours. The coupling constant of effective 5 flavour QCD, 
af^ and the coupling constant of full 6 flavour QCD, af^ both obey the renormalization 
group equation ( with nj = 5 for af^ and rif = 6 for a^ ) 
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| and 

C A -- 3 are the Casimir operators of the fundamental and adjoint representation of the 
colour group SU(3), T F = | is the trace normalization of the fundamental representation. 
The solution of eq. (0) in the next-next-to- leading order has the standard form 



The three loop QCD beta function in the MS'-scheme was calculated in W7\. C F 
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where it is understood that the scale A^^- also depends on the number of active flavours. 



The (top quark mass dependent) relation between a^ and a^ is called the decoupling 
relation and will be discussed later in more detail but for completeness we give the NNL 
order expression 
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where \x is the renormalization scale and m t (fi) is the top quark mass in the MS scheme. 
Please note that the term ^T F C F that we found is slightly different from the one in ref. 

|18| . Substituting expression (H) for a^ and af^ one can find the connection between A 

and A-jro y i a Wt(/i)- 

mt([x) obeys the renormalization group equation 
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and ( is the Riemann zeta-function. The three loop quark mass anomalous dimension was 
calculated in ref. ITS 



Let us quote the existing results for the axial vector contribution. The axial vector 
non-singlet part can be reduced to the vector case by using the effective anticommutation 
property of the 75 matrix in the prescription that we use. More strictly, it can be done 
only in the limit of massless light quarks. Thus the non-singlet axial part coincides with 
the vector non-singlet part up to a change of the weak coupling constants and reads (in 
the leading order of the large top quark mass expansion) 
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Let us turn to the singlet axial vector part. In the Standard Model quarks in a weak 
doublet couple with opposite sign to the Z° boson in the axial vector part of the neutral 
current. That is why the contributions from light doublets add up to zero in the massless 
limit for axial vector singlet diagrams. The only non-zero contribution comes from the 
top-bottom doublet due to the large mass difference between top and bottom quarks. 

The axial vector singlet part in the leading order of the large top quark mass expansion 
has recently been calculated in refs. || [K|. The result in the effective theory with 5 active 
massless quark flavours is 
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where we separated the two weak coupling structures (which was not present in refs. [P, 10]) 
and used the notation g h ^ = g\ and g^ p = g\. Here and below m t = m t (m z ) is the MS 
top mass at the scale mz- One may relate it to the pole mass through the expression 
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which is known in the NNL approximation [^1 or relate it to m t {m t ) through the expression 
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This would correspondingly modify the coefficients of the a s term in ( JT~3|) but we prefer to 
use the MS top quark mass, m 4 (/i) (at /i = m z ) which is the original mass from the QCD 
Lagrangian. 



In the present paper we present the power suppressed top quark mass corrections for 
both the vector and axial vector contributions. The Feynman diagrams that we have to cal- 
culate to obtain the power suppressed top quark mass corrections are given in figures 1,2,3. 



Figure 1. Massive diagrams (= with top quark loops) contributing to the vector 
non-singlet part, T h ' ad and axial non-singlet part, T had . The symbol (8> is used 
to indicate an external vertex of the neutral weak vector current for r, ' ■ and 
an axial vector vertex for T had . It is understood that for each diagram at least 
one fermion loop has to be a massive top quark loop and a loop that contains the 
external vertices is always a massless quark loop. The massless diagrams that were 
already calculated in ref. Q are not considered here. 



Figure 2. Massive diagrams (= with top quark loops) contributing to the vector 
singlet part, T had . The symbol (8> is used to indicate a vector current vertex. It is 
understood that for each diagram at least one fermion loop has to be a massive top 
quark loop. The massless diagrams (of the same topologies) were already calculated 
in ref. 101 and are not considered here. 



Figure 3. Massive and massless diagrams contributing to the axial vector singlet 
part, T had . The symbol (g) is used to indicate an axial vector current vertex. For 
each fermion loop both massless quarks and massive top quarks are considered. 



3 The calculation of the massless diagrams 

In this section we will treat the calculation of the four-loop massless diagrams from figure 3 
that contribute to T h ' ad . We will illustrate the techniques by considering the most difficult 
diagram (the first one on the second line in figure 3). 

Since we are only interested in the structure function Il^g 2 ) (see equations (|]||)), we 
contract the diagrams with the projector (g^ — q fJ- q u /q 2 ) which reduces the diagram to a 
scalar integral. In figure 4 we present diagrammatically the renormalization of ultraviolet 
divergences of this diagram. 



Figure 4. The ultraviolet renormalization of the four-loop diagram. The ultraviolet 
counterterms are presented between round brackets. 

At present, we do not have a technique which would allow a direct calculation of this 
4-loop diagram. But we can use the fact that we need to know only poles in e for this 
diagram, since only these pole terms generate terms containing ln(Q 2 //i 2 ) which produce 
non-zero imaginary parts (logarithms come from i(fe) E = ~ + ln(^2-) + 0(e) and give 
imaginary parts through ln(— s — ie) = ln(s) — iir, where s = —Q 2 = q 2 = m 2 z ). One 
can see from figure 4 that it is sufficient to calculate all renormalization terms [2 nd - 4 th 
terms in figure 4), including the 4-loop counterterm, in order to restore the pole terms for 
the diagram itself (because the sum is finite). All renormalization terms except the 4-loop 
counterterms can be directly calculated with the help of the package MINCER [ET1] written 
for the symbolic manipulation program FORM |22f . This package calculates analytically 



3-loop massless propagator diagrams using the integration by part algorithm of ref. |23[ 
for dimensionally regularized diagrams. 

In this way we have reduced the problem of calculating the imaginary part of the 4-loop 
diagram of figure 4 to the problem of calculating its 4-loop ultraviolet counterterm. The last 
problem can be reduced to the calculation of 3-loop propagator type diagrams (which are 
calculable with the package MINCER) using the infrared rearrangement method [f24], p5[ . 
This method relies heavily on the fact that in the MS'-scheme ultraviolet counterterms 



are polynomials in momenta and masses [pq] , i.e. do not contain logarithms or inverse 



powers of momenta and masses. In our case the 4-loop counterterm has dimension two 
which means that it is simply proportional to Q 2 . If we take the d'Alembertian in Q of 
this counterterm we get a dimensionless quantity which is just a Laurent series in e. To 
obtain this Laurent series we need to calculate counterterms of dimensionless diagrams. 
These diagrams are produced after applying the d'Alembertian in Q to the 4-loop diagram 
of figure 4. In fact, the application of the d'Alembertian produces several dimensionless 



4-loop diagrams after differentiation of the lines of the original diagram; some of these 
diagrams are presented in figure 5. 



Figure 5. Diagrammatic representation of applying the d'Alembertian in Q to the 
4-loop counterterm. A prime on a line denotes the differentiation of this line in its 
momentum. The round brackets denote the ultraviolet counterterms of the diagrams 
inside the brackets. 
Since counterterms of dimensionless diagrams do not depend on the external momen- 
tum Q, we can change the route of this momentum through these diagrams as we wish, 
e.g. as it is chosen in the r.h.s. of figure 5. The choice of a new route for the external 
momentum can generate infrared poles (even though the original diagram did not have 
infrared divergences) which then essentially complicates the extraction of the ultraviolet 
poles. For example, nullifying the momentum Q will nullify the whole diagram because of 
new infrared divergences. The chosen momentum route in the r.h.s. of figure 5 reduces the 
calculation of the 4-loop diagrams in figure 5 to the calculation of simpler 4-loop topolo- 
gies. These topologies have the form of a 3-loop propagator type subdiagram (which can 
be done with the package MINCER) inserted in a trivial one-loop topology. 

The essential complication now is that the second diagram in the r.h.s. of figure 5 
contains infrared divergences. For this diagram it is impossible to choose the route of the 
external momentum in a way that avoids infrared divergences and simultaneously reduces 
the calculation of the corresponding diagram to a 3-loop propagator insertion. In this 
case we need to apply the technique of the R* operation |27| which allows to calculate 



ultraviolet counterterms of dimensionally regularized diagrams even in the presence of 
infrared singularities. It is interesting to note that the diagram of figure 4 was the only 
diagram that needed the application of the R* operation in the calculation of T h ^ d . 
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4 The expansion of massive diagrams 

We will now treat the large mass expansion of the individual massive diagrams. We calcu- 
late all integrals in Euclidean momentum space. The general theory of Euclidean asymp- 



totic expansions was developed in |28], p9| . For practical purpose we use the techniques 
developed in f3(J, |3l] . 



Let's go through some simple ideas (see e.g. P^| ) that generalise to a recipe of expand- 
ing individual dimensionally regularized diagrams. This recipe can then also be used to 
expand MS renormalized diagrams since it can be applied to each term of the renormal- 
ized expression for a given diagram (i.e. the expression after application of the ultraviolet 
R-operation to this diagram). 

A simple scalar diagram containing both massless lines and massive lines is shown in 
the l.h.s. of figure 6a. 



Figure 6a. Thick lines denote massive scalar propagators, thin lines denote massless 
scalar propagators. The symbol x indicates the insertion of the small momentum 
expansion of the massive triangle. 

We are going to expand this diagram in a large mass which is equivalent to the expansion 
in a small (in comparison to the mass) external momentum, Q. Please note that we can not 
simply expand the integrand of the corresponding Feynman integral as a Taylor series in Q, 
because putting Q = generates infrared divergences when one integrates over momentum 
k. Let us therefore first consider the expansion of the massive one-loop subgraph as a Taylor 
series in its external momenta, k and Q. This expansion can be obtained by a simple Taylor 
expansion of its integrand in k and Q because it does not generate infrared divergences 
when one integrates over 1. More generally, the expansion of diagrams with only massive 
propagators in terms of (small) external momenta can be safely done by making Taylor 
expansions in these momenta in the integrands. The Taylor expansions of integrands are 
generated by simple expansions of propagators in small momenta Q 

1 _ 1 f (-2P.Q-QJ ( 1 \ 

(P + Q) 2 + M 2 P 2 + M 2 ^ (P 2 + M 2 ) i \(P 2 + M 2 f N + 2 ) ) 

where N is the desired order in the expansion. We then add and subtract the expansion of 
the massive one-loop subgraph as is indicated in the r.h.s. of figure 6a. 

The first term in the r.h.s. of figure 6a is a massless one-loop integral with the expanded 
massive subintegral inserted in the integrand and it can be evaluated without difficulty, 
keeping Q finite. 
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The second term in the r.h.s. of figure 6a (the combination in the square brackets) has a 
vanishing contribution from the integration region where k is small because the behaviour 
of the massive subintegral for small k is subtracted off until the necessary order (that is 
determined by the depth of the expansion). One may then Taylor expand the integrand 
of the second term around Q = to obtain an integral without an external momentum 
(a tadpole integral) that can also be evaluated without difficulty It should be noted that 
massless tadpole diagrams vanish in dimensional regularization which means that after the 
Taylor expansions in the integrand of the second term in Q, only the two loop massive 
tadpole contribution survives. 

Finally we get the large mass expansion as it is presented in figure 6b where we used the 
notation that a box around a (sub)graph indicates that the integrand of this (sub)graph 
is Taylor expanded up to the desired order in its external momenta. This will be the 
standard notation in the following. Please note that the nullified term in figure 6b is a 
massless tadpole. Although the expansion procedure produces new ultraviolet and infrared 
divergences in separate terms of the r.h.s. in figure 6b, these terms are well defined for 
non-zero e and these new divergences cancel in the sum. 



Figure 6b. The diagrammatic large mass expansion of the diagram 6a. 

Let us consider another diagram presented in figure 7a. The same reasoning as for 
diagram 6a holds but now we should first add and subtract the small momentum expansion 
of the massive propagator as it is shown in the r.h.s. of figure 7a. 



Figure 7a. A scalar triangle diagram with one massive line 

After this we can expand the expression in the square brackets in its external momenta. 
Using again the property that massless tadpoles vanish in dimensional regularization we 
are left with the diagrammatic expansion in figure 7b. Please note that the nullified term 
in figure 7b is a massless tadpole. 
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Figure 7b. The diagrammatic large mass expansion of the diagram 7a. 

The triangle diagram of figure 7 can be part of a larger diagram such as the scalar 
diagram presented in figure 8. For figure 8, the same reasoning as for the diagrams in 
figures 6 and 7 holds and we are left with the indicated diagrammatic expansion. 



Figure 8 The diagrammatic large mass expansion of a 2-loop diagram. 

The large mass expansion of two loop propagator type diagrams is also treated in the 
recent article ref. |33fl . It interesting to note that an analogous reasoning (i.e. focusing on 
infrared regions) can be applied to understand the related problem of small mass expansions 



for which a recipe is given in [BO, 34 



Let us now formulate the recipe for the large mass expansion of diagrams. A line with 
a large mass, M, will be called a heavy line. 

First we have to find all asymptotically irreducible (sub)graphs. An asymptotically ir- 
reducible (sub)graph, gAi, is a connected subgraph which contains at least one heavy line 
(and all heavy lines connected to this one via heavy lines) and which can not be made 
disconnected by cutting a non-heavy line. We have to expand each asymptotically irre- 
ducible (sub)graph as a Taylor series in terms of its external momenta. We graphically 
indicate this by drawing a box around the asymptotically irreducible (sub)graph. Then 
the large mass expansion of the whole Feynman diagram is the sum over all combinations 
of non-overlapping boxes that can be drawn in this diagram so that all heavy lines are in 
boxes. Note that a box does never cut a heavy line. 



Or in a symbolic form: 

G= E % f ) *T (A0 W + O(^), (14) 

{9Al} 

where G is the Feynman graph to be expanded. The sum goes over all sets, {gAi}, of non- 
overlaping asymptotically irreducible (sub)graphs comprising all heavy lines. T^ {gAi} 
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denotes the Taylor expansion of (sub)graphs from the set {gAi} in their external momenta 
until the necessary order. G/{gAi} is the graph obtained from G by shrinking the subgraphs 
from {qai} into points. 

As a practical example we will now treat the diagrams that contribute to the axial 
vector part of the Z boson decay rate. At the 3-loop level the only Feynman diagrams 
that contribute are so called 'double triangle diagrams' (the first topology in figure 3) 
with the triangles formed by bottom and top fermion loops. This correction was originally 
calculated in || and confirmed in J7|. The diagram with two massive triangles is zero 
since the only physical cut is through two gluon propagators and Z° decay in two gluons is 
kinematically forbidden (Landau- Yang theorem). We are therefore left with two diagrams 
to be calculated, one massive and one massless diagram. The massless diagram can be 
directly calculated with the package MINCER. The massive diagram is calculated with the 
above recipe. 

A diagrammatic representation of the ultraviolet R operation followed by the asymp- 
totic expansion procedure, applied to the massive double triangle diagram is given in 
figure 9. 



Figure 9. Thick lines indicate top quark propagators, thin lines indicate the mass- 
less quark propagators and spiral lines - gluon propagators. Ultraviolet counterterms 
are indicated between round brackets. The asymptotically irreducible (sub)graphs 
are surrounded by boxes with the corresponding tadpole topologies indicated below. 
After the Taylor expansions the subgraphs in boxes are reduced to massive vacuum inte- 
grals and the resulting master topologies are indicated under the boxes. When the boxed 
subgraphs are integrated out we are left with massless diagrams that can be calculated 
with the package MINCER. We have written efficient FORM procedures to perform the 
necessary massive vacuum integrals. The procedures include 3-loop topologies of the type 
Benz and non-planar and use recursions based on the recursion scheme of |)35|| . They are 
essential for the large mass expansion of 4-loop diagrams. We should emphasise that only 
with very efficient massive vacuum procedures, that can deal with tensor numerators, the 
4-loop calculation is feasible. 
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Making the Taylor expansions deep enough, we find the results for the contributions in 
figure 9. Adding the massless double triangle diagram and taking the imaginary part by 
applying ln(— s — ie) = ln(s) — in (remember that s = m 2 z is the squared 4-momentum of 
the Z boson) we reproduced the large mass expansion of the known result || 

r£Z = ?# (-Y U + 4- 2 + d 2 A + dl^) (is) 



8v/2' 



7T 



with <f 2 = T 2 D [-g + |ln(^)], d\ = T 2 D (£), d\ = T 2 D (^), d\ = T 2 D (- 



52 



,165375 

D = n 2 — 1 is the number of generators of the colour group SU(n) (D = 8 for QCD). 

We will now show that the same method can be used to compute the a^ correction to 
T had . In contrast to the previous 3-loop case where we could determine both the real and 
imaginary parts of a diagram, we are now able to compute only the necessary imaginary 
part of diagrams with the present available techniques (we don't have analytical results for 
general 4-loop massive tadpoles). In order to find the imaginary part of a diagram we need 
to calculate only terms that contain logarithms of s (since we take the imaginary part by 
applying ln(— s — ie) = ln(s) —in). In dimensional regularization every massless propagator 

2 2 

type integral receives a factor (— ) e = 1 + eln(^-) + 0(e 2 ). In contrast, vacuum massive 

2 

integrals produce factors (-^) £ (m t is the top quark mass) and do not give logarithms of 
s. This means that contributions without massless propagator type integrals do not have 
to be considered for our purpose. The expansion of one renormalized 4-loop diagram that 
contributes to T h ^ d is presented in figure 10. 



Figure 10. The large mass expansion of a renormalized 4-loop diagram 

In figure 10 we have only presented terms which contribute to the imaginary part of 
the diagram. All 4-loop massive diagrams contributing to r^ give after the large mass 
expansion topologies that can be calculated with MINCER (the massless parts) and with 
the tadpole procedures (the vacuum massive parts). 
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5 The results for T 



had 



The results of the a 3 , approximation for T h ' ad , T h ^ d , r h ' ad and T had in the order 1/m® of 
the large top quark mass expansion are presented below. Note that these results are in 
6-flavour QCD, Nf = 6 and the decoupling relation (|9|) has not been applied yet. 
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where n is the parameter of the colour group SU(n) (n = 3 for QCD). Note that the 
coefficient b\ agrees with ref. |Rj| and 6 2 , 63 agree with the expansion of the exact result 
ref. 1371. 
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For completeness we also present separately the bottom-bottom and top-top contribu- 
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weak coupling constants: the bottom-bottom, top-top and top-bottom contributions. 
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Note that except for the massless diagrams there is only one masssive diagram that con- 
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tributes to T had ' (the 3 topology in the 3 line of figure 3) 
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Note also that although separate diagrams that contribute to T had ' are generally non- 
zero, they cancel in the sum except for the one diagram (the 3 rd topology in the 3 rd line 
of figure 3) that has a colour factor proportional to T F D(Nf — 1). This fact can be un- 



derstood through the operator product expansion technique ref. |^TJ. In addition note 
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that the imaginary part of the diagram with three top quark loops (the 3 rd topology in 
the 3 rd line of figure 3) vanishes according to the Landau- Yang theorem (which we checked 
explicitly). 

The result for the vector singlet part reads 
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where <i a6c are the symmetrical structure constants of the SU(n) colour group, d abc d a b c = 
40/3 for QCD. Note that the constant c 3 comes from massless diagrams only (although 
separate massive diagrams are non-zero in the leading order of the large top mass expan- 
sion, they add up to zero) as it is required by the decoupling mechanism for the vector 
part. It is interesting to mention that the diagrams with two massive top quark loops are 
non-zero separately but they add up to zero in all orders of the large mass expansion. 
Explicit checks show that the coefficients of the logarithms in eqs. (16) and (|Tfl) are 



in agreement with the required renormalization group invariance of the physical quantity 
which in the a 3 , approximation reads 

V 2 -^- 2 T(^,a s (^m t ^)) = 0(at). (21) 

Of course the true physical quantity is T had = T had + T had + T had + T had . But from a 
theoretical point of view each of the four separate parts is renormalized independently and 
is therefore renormalization group invariant by itself. 

The results that are presented in this section were obtained in an arbitrary covariant 
gauge for the gluon fields i.e. keeping the gauge parameter as a free parameter in the calcu- 
lations. The explicit cancellation of the gauge dependence in the physical quantities gives 
a good check of the results. Individual diagrams contain £2, and £4 but these contributions 
add up to zero in the total results. 
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6 Final results (in a decoupled form) 

It is known that the Appelquist-Carazzone theorem |3£| about the decoupling of a heavy 
particle in quantum field theory does not work in its naive form for the MS renormalization 
scheme and one should make an extra shift in the coupling constant (see eq.(|9]) ) to make 
the decoupling explicit (i.e. to kill the large-mass logarithms) |3^, [l£|. However in the 
presence of an axial vector current the decoupling does not work (even after the shift in 
the coupling constant) due to the presence of axial anomaly type diagrams as is the case 
for the axial vector singlet contribution to T^ad- 

Since we have explicitly calculated the top quark mass terms in the order a^ for T had , we 
can derive the decoupling relation for the QCD coupling constant in the next-to-next-to- 
leading (NNL) order. Because the vector contribution to Thad should obey the decoupling 
mechanism, the use of the decoupling relation should convert the new 6 flavour result for 
T had (see eq.( [L6|) ) to the previously known result in effective massless 5 flavour QCD (see 
eq.(||)). One can see that the NNL order decoupling relation obtained in this way (see 
eq.(|9])) slightly differs from the one known in literature [HI 



In order to settle this discrepancy, we performed an independent calculation. We have 
calculated the 3-loop massless quark propagator with a zero momentum ipip operator in- 
sertion (where ip is the quark field) 

G W ^(Q 2 ) = J e^dxdy < o|T {f a (x)My)My)M®V s } 10 > (22) 

We will derive the decoupling relations from this (gauge dependent) Green function. One 
can also use the normal quark propagator for this purpose but one has to evaluate it at 
4-loop level to derive the decoupling relation for the coupling constant in the NNL order. 
This is because at the one loop level the quark propagator is proportional to the gauge 
parameter. 

For a physical quantity the decoupling mechanism consists of a shift in the coupling 
constants (and in the light masses, if present). The decoupling mechanism for the Green 
function is more complicated than that for a physical quantity because the renormalization 
of the Green function involves an overall renormalization constant, G ren = ZGsare, and 
one should also obtain the decoupling relation for Z. To avoid this small complication, we 
prefer to take the quantity 



d\n(Q 

that imitates a physical quantity in the sense that it has no overall renormalization constant 
Z. The result of our calculation in the leading order of the large top quark mass expansion 
reads 



21 



dln(Q 2 



In (G w 



ip[ipil>]ip,ren 



{Q 2 



a 



(6)' 



7T 



0,(1 " 3« 



+ 



o 



(6)' 



7T 



N f T F C F 



5 1 Q 2 , 
9 + 3 ln( ^ 



H-TpCf 



3 ra| 



+CaC f 



295 11 ,Q* , 



1 
6 



i, ( Q 2 - 



+ ^ 



1 

48 



+ 



a 



(6)' 



7T 



N f T F C F 



2 , ^%l (^V 10 i rA 



9 



m? 



+^C^J^ + ^ln(^)ln(% 
19 18 mf ;U 2 



2 



+£ 



89 



1728 24 



lln(4)ln(^) + 

mi jj. 



385 /*■ 

216 nl m t 2J 

13 , / A ! i 2 

144 mf 48 



2 ■ 



+T F Cf 



13 1 , ,/i 2 . 
"48 + 12 ln( ^ 



1 
9 



F T^ln 2 (^) 



o 



(6)' 



7T 



-M— 



24 



massless 
contrib. 



+ Ci 



(23) 



up to an overall normalization. £ is the gauge parameter that appears in the gluon propa- 
gator as 2 l +ifL {—g^ u + (1 — £) g 2+i e )- For the af term only contributions that contain one or 
more top quark loops are presented. These massive contributions of the order a^ should 
disappear after the application of the decoupling relations for the coupling constant and for 
the gauge parameter £. From this requirement we find the decoupling relations in the NNL 
order for the coupling constant presented in eq. @ and for the gauge parameter presented 
below. 
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Please note that the term -^T F C F that we found in eq.(|4|) for the gauge parameter (as 
well as the analogous term for the coupling constant) is slightly different from the one in 
ref. fL8"f . Note that the decoupling relation for the coupling constant, eq.@, is derived by 
us in two independent ways: from the 4-loop calculation of Thad and the 3-loop calculation 
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We will now give the results for Thad in effective QCD with 5 active massless flavours 
at the renormalization scale \x = mz- These results are obtained by substitution of the 
decoupling relation, eq.(H), into the results for 6 flavour QCD of eqs.(|T6|), (|TTD and (|20|). 
We will use the notation x 
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6 2 = 1.4092 

+ [0.065185 - 0.014815 ln(x)] x 

+ [-0.0012311 + 0.00039683 ln(x)] a; 2 

+ [0.000061327 - 0.000023516 ln(x)] x 3 + 0(x 4 ) 
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+ 
+ 
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-0.17374 + 0.21242 ln(x) -0.037243 In 2 (x)\ x 
-0.0075218 - 0.00058859 ln(x) + 0.00038305 ln 2 (x)] x 2 
0.00050411 - 0.00012099 ln(ar) + 0.000031419 ln 2 (x)l x 3 + 0(x 
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.r 2 



+ [-0.0031322 + 0.012117 ln(x) - 0.0011905 In 2 (x) 

+ -0.00088827 + 0.00047262 ln(x) -0.00017637 ln 2 (x)] x 3 + 0(x 4 ). 

We want to stress that the results show an excellent convergence of the large top quark 
mass expansion which can be seen from the fast decrease of the coefficients with the order 
of the expansion parameter x = ^f . 

Note that the massive logarithms ln(^f) are present in the leading order of the large 

mass expansion in T had only (the violation of decoupling of the top quark). These log- 
arithms can be summed up using the operator product expansion technique to produce 
a result that is finite in the limit of an infinitely large top quark mass as it was done in 
ref. for the order a 2 . However, for realistic values of m t the above expressions can be 
trusted and are quite stable with respect to a change in the renormalization parameter /i 
around the natural scale for this process /i = mz, as for example can be seen for T had in 
figure 11. 



A Q 

Figure 11. The u dependence of T had around u = mz for m t (mz) = 140 GeV and 

(5) 

^-QCD = 0-2 GeV. The dotted line indicates the 3-loop result, the solid line indicates 
the result up to (and including) 4-loops. 

One can use T had from the present paper to obtain the hadronic decay width r(W /± — > 
hadrons) (by a standard change of the weak coupling constants) and to obtain the large 
mass expansion in the a^ order for r(r~ — ► v T + hadrons) (as we do it in the next section). 
From the results of this paper one can also straightforwardly obtain the a^ approximation 
to the total cross section of electron-positron annihilation a tot (e + e~ — > 7, Z° — ► hadrons) 
in the necessary energy range, e.g. below, or above the Z° peak. 

24 



7 The large charm quark mass expansion for r lepton 
decay 

Another prominent process (beside Z boson decay) to extract the value of a s from ex- 
periment is semihadronic r lepton decay. Although the scale of this process is relatively 
low, non-perturbative corrections turn out to be small and perturbative QCD can be used 
ICJ-fJ^l to calculate the r lepton decay rate, or the ratio 
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v T e v f 



The ratio R T is expressed through the imaginary part of the W boson current correlator 
as (modulo a contribution that vanishes for massless active flavours) 
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where the integration is over the invariant mass of the hadrons, s, and we use the normal- 



izations of ref. j|!|. II! is the transverse part of the W boson current correlator analo- 
gous to the Z boson current correlator of eq. @ (in eq. @ the neutral weak quark current 
should be replaced by the charged weak quark current J M = -2^ J2i,j uil^ (1 — 7 5 ) Vijdj 
with Ui = (u, c, t) and dj = (d, s,b)). 

We work in the approximation of massless u, d, s quarks and a heavy c quark within 
perturbative QCD omitting non-QCD corrections. The a^ approximation in the leading 
order of a large mass expansion was calculated in ref. 0. The first charm quark mass 
suppressed term of the order a 2 s is obtained in ref. p6| . We will obtain the large c quark 

2 

mass expansion of R T in the parameter ^j within effective QCD with 3 massless flavours. 
This is the correct expansion parameter since the mass of the r lepton is below the threshold 
for the production of charmed hadrons. Therefore a charm quark appears only in internal 
fermion loops, the effective expansion parameter appears to be m^/(4m^) and a large c 
quark mass expansion is justified (although the charm quark mass m c is smaller that the 
tau lepton mass m T ). The Feynman diagrams that contribute to ImUi for the W boson 
are of the non-singlet type only (see figure 1) and were already calculated for the case of 
Z boson decay (see eq.(|16|)). After performing the integration in eq.fl2"9|) with the results 
of eq. fll6l) we have to apply the decoupling relation to go to effective QCD with 3 active 
flavours. Putting \i = m T we obtain 
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where ai 3 -* = a^\ 



m T 



m c = m c {m T ) is the MS charm quark mass and n = 3 is the 



number of quark colours. We neglected terms that are suppressed by the bottom and top 
quark masses (which in principle are present after the decoupling). Substitution of the 
QCD colour factors and rif—3 gives 
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J 
,2- 



1, r° = 5.2023, 

0.023778 - 0.0044444 In (^|) , 

-0.00020118 + 0.000052910 In (^|) 

0.0000053203 - 0.0000016797 In (^f ) , 

26.3659, 

-0.057156 + 0.079881 In (^4) - 0.012654 In 2 f =4 

-0.00099668 - 0.00016858 In (2j) + 0.0000687096 In 2 (=4 

0.000036522 -0.0000089407 In (=4) + 0.00000168435 In 2 (4- 

V !ri c / V rn c 



Note that the coefficient r\ agrees with ref. [36]. 

2 / \ 2 

Although the expansion parameter ^j p» ( xz±a f 3 ) ^ s slightly larger than 1, the fast 
decrease of the coefficients ensures a good convergence of the large charm quark mass 
expansion. 

We conclude that the large mass expansion converges fast for both Z boson and r lepton 
decays and the obtained a\ approximations can be trusted. 
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